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ABSTRACT 

We explore the kinematic properties of dense continuum clumps in the Perseus molecular 
cloud, derived from our wide-field C^^O J = 3 ^ 2 data across four regions - NGC 1333, 
IC348/HH211, L1448 and L1455. Two distinct populations are examined, identified using 
the automated algorithms CLFIND (85 clumps) and GAUSSCLUMPS (122 clumps) on exist- 
ing SCUBA 850 |Ltm data. These kinematic signatures are compared to the clumps' dust 
continuum properties. We calculate each clump's non-thermal linewidth and virial mass 
from the associated C^^O/ = 3^2 spectrum. The clumps have supersonic linewidths, 
(ontM) = 1.76 ±0.09 (CLFIND population) and 1.71 ±0.05 (with GAUSSCLUMPS). The 
linewidth distributions suggest the C^^O line probes a lower-density 'envelope' rather than a 
dense inner core. Similar linewidth distributions for protostellar and starless clumps implies 
protostars do not have a significant impact on their immediate environment. The proximity to 
an active young stellar cluster seems to affect the linewidths: those in NGC 1333 are greater 
than elsewhere. In IC348 the proximity to the old IR cluster has little influence, with the 
linewidths being the smallest of all. The virial analysis suggests that the clumps are bound 
and close to equipartition, with virial masses similar to the masses derived from the con- 
tinuum emission. In particular, the starless clumps occupy the same parameter space as the 
protostars, suggesting they are true stellar precursors and will go on to form stars. We also 
search for ordered C^^O velocity gradients across the face of each core. Approximately one 
third have significant detections, which we mainly interpret in terms of rotation. However, we 
note a correlation between the directions of the identified gradients and outflows across the 
protostars, indicating we may not have a purely rotational signature. The fitted gradients are in 
the range ^ = 1 to 16 km s~^ pc~^ larger than found in previous work, probably as a result of 
the higher resolution of our data and/or outflow contamination. These gradients, if interpreted 
solely in terms of rotation, suggest that the rotation is not dynamically significant: the ratios of 
clump rotational to gravitational energy are jSrot ^ 0.02. Furthermore, derived specific angular 
momenta are smaller than observed in previous studies, centred around j ^ 10~^ kms~^ pc, 
which indicates we have identified lower levels of rotation, or that theC^^O/ = 3^2 line 
probes conditions significantly denser and/or colder than ^ '^ 10^ cm~^ and T ^ 10 K. 

Key words: submillimetre - stars: formation - stars: evolution - ISM: kinematics and dy- 
namics - ISM: individual: Perseus. 



1 INTRODUCTION 

Stars form inside dense cores deep within molecular clouds. The 
similarity of the core mass function (CMF) to the initial mass func- 
tion of stars (IMF, see |Motte, Andre & Neri| 1998 , Alves, Lombardi 
'& Lada'2007) has sometimes been used (e.g. Enoch et al. 2008 ) as 
evidence against theories of core formation where no correspon- 



E-mail: e.curtis@mrao.cam.ac.uk 



t E-mail: jsriUis^cam.ac.uK 



dence is anticipated, e.g. in the competitive accretion picture ( |Bon-| 
|nell et al.|20"0T]|Bate & BonnelipOOg] ). However, the mapping of a 
given CMF on to the resultant IMF is complicated by many factors 
(see e.g. Hatchell & Fuller 2008, Curtis & Richer 2010), such as 
the cores' varying multiplicity (Good win et al.,2008^ , star-forming 
efficiencies ( Swift & Williams 2008) and/or mass-dependent life- 
times (Clark, Klessen & Bonnell 200 7}. In fact, as s hown by | Swift] 
|& Williams ( 2008 ) and G oodwin & Kouwenhoven| ( |2009| >, diverse 
evolutionary schemes can map the observed CMFs on to the IMF, 



© 2010 RAS 



E, L Curtis and J, S. Richer 



implying the current mass function data are not sufficient to dis- 
criminate between different theoretical models by themselves. 

The kinematics of star-forming cores, probed with spectral- 
line observations of dense molecular tracers, are arguably the best 
discriminator between different models of core formation (e.g. 
[Andre et al.|2007l|Kirk, Johnstone & T afalla 2007). For example, 
cores created by shocks in large-scale flows exhibit large veloc- 
ity gradients and are located at local maxima in the line-of- sight 
velocity dispersi on (e.g.|Ballesteros-Paredes, Klessen & Vazquez-] 
|Semadeni||2003|[Klessen et al.||2005| ). Alternatively, the magneti- 
cally controlled scenario has more quiescent velocity fields (e.g. 
[Nakamura & Li 2005 ), well-matched to observations of isolated 
starless cores with subsonic levels of turbulence (e.g. i Myers|1983l 
[Caselli et al.|2002| ). Such quiescent cores are in opposition to purely 
hydrodynamic models of gravoturbulent fragmentation (see "Mac" 
[Low & Klessen|2004| ) which produce a majority of cores with su- 
personic velocity dispersions. 



1.1 Nomenclature and dump populations 



We adhere to the clu mp hierarchy explained in [Cu rtis & Richerj 
( |20T0| ), which follows [WTiams, Blitz & McKeeJ pOOO). Briefly, 
molecular clouds contain clumps which in turn harbour cores, the 
direct precursors of individual or multiple stars. Clumps with- 
out embedded objects are starless, unless they are gravitationally 
bound when we refer to them as prestellar. 



2 OBSERVATIONAL DATA 

2.1 Dust continuum 

The SCUBA 850 |^m data analysed have a beam size of 14 arcsec 
(0.017 pc at 250 pc, our adopted distance to Perseus, see [Paper T] l 
and are sampled on a 3 arcsec grid with a ty pical rms noise level, 
arms = 35mJybeam"^ JHatchell et al.|2005| ). 



[Offner, Klein & McKee[ ( |2008a| ) and [Krumholz, McKee &[ 
[Klein[ ( [2005| ) maintain that the physical mechanism of star forma- 
tion depends on the much debated presence (or absence) of turbu- 
lent feedback. In one view, molecular clouds are short-lived (on 
timescales of order one dynamical time), non-equilibrium struc- 
tures (e.g. [Elmegreen et al.|2000[[Hartmann 2001 , Dib et al.|2007] ) 
characterized by transient turbulence, which dissipates quickly. In 
this scenario, stars could form by the collapse of discrete cores. 
Conversely, molecular clouds might form slowly and be quasi- 
equilibrium objects (see e.g. Shu, Adams & Lizano'1 987|[McKee[ 
[1^991 [Krumholz & Tan 2007 , Nakamura & Li 2007). This would 
require turbulence to be constantly injected into the clouds, either 
externally (from e.g. supernova blast waves or Hll regions) or inter- 
nally (from e.g. protostellar outflows) and cores could form through 
competitive accretion ( Bonnell et al.[[2001| l. [Offner et aL[ ( |2008b| l 
demonstrate that these two disparate views of star formation give 
rise to distinguishable kinematics in dense stellar precursors and 
thus the study of core kinematics potentially offers a way to probe 
the turbulent state of the surrounding molecular cloud. 

This paper is concerned with the kinematical properties of 
dense clumps in the Perseus molecular cloud (hereafter simply 
Perseus) and compares them to theoretical models. We examine 
the linewidths, virial masses and ordered velocity gradients of two 
populations of clumps (Curtis & Richer 2010) identified across 
SCUBA dust continuum maps ( ,Hatchell et al.,2005^ . These kine- 
matic signatures are derived for each clump from our wide-field 
survey ( Curtis, Richer & Buckle||2010a[ [Paper I I of four clusters 
of star-forming cores towards NGC 1333, IC348/HH211 (simply 
IC348 hereafter), L1448 and L1455 in the 7 = 3 ^ 2 rotational 
lines of CO and its common isotopologues ^^CO and C^^O. In 
the preceding paper of this series ( jCurtis et al.|[2010bl [Paper II| ), 
we explored the evolution of molecular outflows across these re- 
gions. We organize this paper as follows: the remainder of the in- 
troduction explains the naming conventions we have adopted to de- 
scribe star-forming clumps/cores, ^briefly describes the contin- 
uum data and associated clump catalogues alongside our spectral- 
line datasets, from which we derive the clump kinematics. Vari- 
ous clump kinematic signatures are investigated inj3] including 
the clumps' linewidths ( §3.1| l and virial masses ( §3. 2| )7 before we 
search for ordered velocity gradients across the face of the clumps 
in §3.3| which we interpret with a solid-body-rotation model. Fi- 
nally, we summarize this work in El 



2.2 Clump catalogues 

We investigate the kinematics of two populations of dust contin- 
uum clumps, which we identified in the SCUBA 850 ^im data in a 
previous paper ( [Curtis & Ri cher 201 0|, using the two most popular 
automated algorithm s: CLFIND (.Williams, de Geus & Blitz|1994| l 
and GAUSSCLUMPS ( [Stutzki & Gusten[|1990| ). Our motivation in 
that study was to determine which clump properties can be robustly 
measured and which depend on the extraction technique. We con- 
tinue that approach in this work, analysing the kinematic signatures 
for each population in turn. 

In [Curtis & Richer [ ( |20T0] ), we located 85 and 122 clumps us- 
ing CLFIND and GAUSSCLUMPS respectively with peak flux densi- 
ties ^ 4(7rms = 140 mJy beam" ^, across the same areas where we 
also have C^^O/ = 3^2 data |Paper l| . The clumps were iden- 
tified as starless. Class or Class I proto stars using the classifica- 
tions of [Hatchell et al.[ ( [?007| hereafter [HFR07| ) , which are based on 
source SEDs incorporating Spitzer data. We associate a clump with 
a HFR07 source, if the separation between the clump and[HFRQ7 



Curtis 



source peak positions is less than the clump's diameter (see 
|&Richer|2010| ). Given the similarities in the extraction method and 
resulting population, most of the CLFIND clumps are directly on top 
of the [HFR07[ sources and are unambiguously identified. However, 
three GAUSSCLUMPS sources are large enough to encapsulate two 
|HFR07[ sources and their results are correspondingly included in 
the population statistics twice. 



2.3 Spectral-line 

The C^^OJ = 3^2 observations used in this analysis were under- 
taken with HARP (the Heterodyne Array Receiver Project; [Buckle| 
|et al.|2009| ) on the James Clerk Maxwell Telescope and have been 
described in detail previously (Paper I). Striping artefacts resulting 
from systematic differences in calibration between HARP channels 
were removed using a 'flatfield' procedure (fully described in |Paper| 
|l]). The final datasets are sampled on a 3 arcsec grid, distributed us- 
ing a Gaussian gridding kernel with a 9 arcsec full- width half max- 
imum (FWHM), yielding an equivalent beam size of 17.7 arcsec. 
The median rms spectral noise values, averaged across each re- 
gion's map (measured in 0.15 kms~^ channels), are 0.18, 0.15, 
0.15 and 0.14 K for NGC 1333, IC348, L1448 and L1455 respec- 
tively. 



© 2010 RAS, MNRAS 000,[l]{T7] 



The kinematics of SCUBA clumps in Perseus 3 



3 THE KINEMATICS OF THE SCUBA CORES 

In molecular clouds, theC^^O/ = 3^2 line is expected to be opti- 
cally thin and excited in dense regions («crit ~ 10^ cm~^). It should 
therefore trace material in the vicinity of or within the dense ag- 
glomerations surveyed wit h SC UBA. At the peak of every clump 
in our two catalogues (see { 2.2 1, we extracted the C^^O / = 3 ^ 2 
spectrum and fitted a Gaussian profile to it using SPLAlQ if the 
peak of the C^^O line was ^ 3(7sp, where CTsp is the rms spectral 
noise. We did not detect C^^O at the 3(7 level towards 4 (5 per cent) 
CLFIND and 6 (5 per cent) GAUSSCLUMPS positions. Of the detec- 
tions, 13 (15 per cent) CLFIND and 18 (15 per cent) GAUSSCLUMPS 
spectra have double- peaked line profiles. For comparison, |Kirk| 



etal 
^0/ 



(2007| hereafter [KJT07| made pointed N2H+ 7 = 1 ^ and 
1 observations on the Institut de Radio Astronomic 



C^«0/ = 2 

Millimetric (IRAM) 30-m telescope (with beam sizes of 25 and 
llarcsec respectively) towards 157 candidate cores in Perseus - 
44 pointings in our fields - and found 66 (42 per cent) C^^O spec- 
tra required two component Gaussian fits with three in one case. 
Generally, / = 2 ^ 1 data trace lower density («crit ~ 10^ cm~^) 
and/or colder material than the 3 ^ 2 at lower optical depths by a 
factor of 2-3 in the optically thin limit. Therefore the / = 2 ^ 1 
line is unlikely to saturate if the 7 = 3 ^ 2 data are optically thin. 
This suggests there are separate emitting sources along the line- 
of- sight with different temperatures and/or densities which require 
multiple-component fits for one transition and not the other. An al- 
ternative is simply that at our velocity res olution (0.15 km s~^ com- 

s separate compo- 



KJT07 



pared to KJTOT's 0.05 km s"^) many of[] 

nents are blended into one, which might partly explain our larger 

line widths (see §3.1| l. 



3.1 Linewidths 

The linewidths in molecular clouds are substantially larger than 
the sound speed of their constituent gas (e.g. jZuckerman & Evans] 
|1974[[So"lomon et al.|1987| l. Such supersonic linewidths are thought 
to arise from chaotic and/or turbulent motions of unknown funda- 
mental origin (see' Elmegreen & Sc alo 2004 and references therein). 
The velocity dispersion in a region scales with its size (e.g. |Larson| 
|1981| ) akin to the Kolmogorov law for subsonic incompressible tur- 
bulence, although turbulence in the interstellar medium is actually 
highly compressible and often supersonic causing shocks. How- 
ever, dense star-forming cores harboured in cloud interiors pos- 
sess near-thermal linewidths almost devoid of the turbulent move- 



ments that cause the broadening in their outer envelopes (e.g. Ben- 



son & Myers 1989"Bai Tanco & Goodman| 1998 [[Goodman et a" 
1998, Pineda et al.,,2010t . The C^^O J = 3 ^ 2 transition is ex 



pected to probe a region intermediate in scale between the dense 
core and its environment, although it is not the tracer of choice 
for young prestellar cores. In the cold (< 10 K), high-density (^ 
10^ cm~^) centres of such objects, the conditions may result in 
depletion, where many molecules (including CO and its isotopo- 
logues) freeze-out on to dust grains. Such zones are better probed 
by molecules such as N2H+and H2D+, but in any case the CO 
J = 3 ^ 2 transitions should probe higher densities than lower / 
ones. 

We take a core's total one-dimensional velocity dispersion, 
(7(-i8o to ^^ the quadrature sum of its thermal and non-thermal dis- 
persions. Or and Ont respectively: 



Table 1. Summary of the turbulent fractions, /turb. derived from Gaussian 
fits to the C^^O / = 3 ^ 2 line at the peak of every SCUBA clump identified 
with CLFIND ( [Curtis & R icher 2010). The results are broken into either 
different regions in Perseus or clump classifications (see p. 2) . The quoted 
errors (cr) are errors on the mean (a/y/N) not sample deviations. 



Population Number /tui-b ct/ 



All 

NGC1333 

IC348 

L1448 

L1455 



81 1.76 0.09 

37 2.05 0.12 

23 1.36 0.09 

14 1.8 0.2 

7 1.5 0.4 



Starless 


17 


1.8 


0.2 


Protostars^ 


33 


1.89 


0.12 


Class 


19 


2.04 


0.16 


Class I 


14 


1.7 


0.2 



^ The Class and I populations combined. 



^Ci80 



2 2 



To find Ont (plotted in Fig.[T]), we calculate Or: 



kT 



HQiSQmp 



(1) 



(2) 



where rup is the mass of a proton and m^isQ the relative molecular 
mass of C^^O (m^iso ^ ^0)- The clump temperature, T, is taken 
from the kinetic temperatures of ammonia gas ( [Rosolowsky et al.[ 
[2008[ ), where available, or 10 K and 15 K for identified starless and 
protostellar cores respectively where not. For those clumps without 
a HFR07 association we assume a temperature of 12K. We calcu- 
late the turbulent fraction, plotted in Figs. [2] to [5] /t^rb = c^nt/^^s. 
where c^ = ^ykT/mmp is the thermal sound speed in the bulk 
of the gas of mean molecular mass m = 2.33mp (assuming 1 
He for every 5 H2) and 7 the adiabatic index (~ 7/5 for a di- 
atomic molecule). For T = \0K, c^ = 0.22 km s~^ and at 15 K 
Cs = 0.27 km s~^. Summaries of /t^rb for various populations are 
listed in Tabs. [T]and[2] 

A majority of the linewidths are supersonic with means 
(/turb) = 1.76±0.0^and 1.71 ±0.05 for the CLFIND and GAUSS- 
CLUMPS populations respectively. The level of non-thermal broad- 
ening should take into account motions from both: (a) the source 
itself i.e. infall, rotation and outflows and (b) the source's natal 
environment. Young stars can inject significant energy into their 
surroundings, increasing the non-thermal motions within nearby 
star-forming cores. For example [Caselli & Myers[ ( [r995[ ) found an 
inverse relationship between a core's NH3 linewidth and the dis- 
tance to the nearest young stellar cluster in Orion B. There are three 
such young clusters in Perseus (see [Hatchell et al.|[ 2005): IC348, 
NGC 1333 and the Per OB2 association. Therefore, of our fields, 
NGC 1333 and IC348 are very near to young clusters, while L1448 
and LI 45 5 are further away (LI 45 5 being the furthest). 

The variations in /t^j-b region-to-region (see Fig.[2|, which we 
might attribute to each region's distance from a young stellar clus- 
ter, are slightly different for the two populations. For the CLFIND 
clumps (see Tab.[TJ, NGC 1333 and L1448 have statistically simi- 
lar fractions: 2.05 ± 0. 12 and 1 .8 ± 0.2 respectively. We ignore the 



^ Part of the Starlink software collection, see |http : //star link . jach7[ 
[hawaii.edu. [ 



^ Errors quoted on averages throughout this paper are errors on the mean 
(a/^N) and not simple standard deviations (a). 
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Figure 1. Spatial distribution of C^^O non-thermal linewidth a cross NGC 1333 (upper left), IC348 (upper right), L1448 (lower right) and L1455 (lower left). 
The grey-scale is SCUBA 850 [xm flux density in mJy beam-^ JHatchell et ari2 005 ) with co ntours at 100, 2 00, 400, 800, 1600 and 3200 mJy beam" i. Circles 
are positioned at the peak of every SCUBA clump in our catalogue identified using CLFIND jCurtis & Richer|2010^ where there is a 3(T C^^O detection. The 
circles' size and colour reflect the linewidth as shown in the key. 



Table 2. Summary of the turbulent fractions, /turb, for the GAUSSCLUMPS 
clump population as in Tab.Q] 



Population 


Number 


/turb 


^/ 


All 


119 


1.71 


0.05 


NGC 1333 


64 


1.91 


0.07 


IC348 


25 


1.44 


0.07 


L1448 


23 


1.52 


0.13 


L1455 


7 


1.5 


0.3 


Starless 


24 


1.61 


0.11 


Protostars 


40 


1.80 


0.10 


Class 


27 


1.73 


0.11 


Class I 


13 


2.0 


0.2 



L1455 populations in these comparisons as they are so few in num- 
ber. IC348, on the other hand, has much smaller linewidths over- 
all: (/turb) = 1.36 ±0.09. The GAUSSCLUMPS population (see Tab. 



|2| has clumps in NGC 1333 with much larger turbulent fractions 
than the similar /t^rb means in IC348 and L1448: 1 .91 ± 0.07 com- 
pared to 1.44 ±0.07 and 1.52 ±0.13. Therefore, for both popu- 
lations, the IC348 linewidths are smaller than in NGC 1333, even 
though both regions are close to stellar clusters. This is possibly 
(as noted in Paper I) because the IC348 IR cluster is old and no 
longer actively forming stars, so it affects its environs much less 
than the NGC 1333 cluster. There is a wide range of clump /t^j-b in 
L1448, which could be a result of a significant contribution to the 
C^^O linewidth from molecular outflows. Many of the protostellar 
clumps in L1448 drive particularly strong flows ( ,Paper II ). 

On dividing the clumps into protostellar and starless subsets 
(see Fig.|3]), there appears to little difference between the distribu- 
tions of turbulent fraction for both source types and clump pop- 
ulations. At the peak positions of the CLFIND objects, we find 
(/turb) = 1 .8 ± 0.2 for starless clumps compared to 1 .89 ± 0. 12 for 
protostars. Similarly, (/^rb) = 1.61 ±0.11 and 1.80±0.10 for star- 
less and protostellar GAUSSCLUMPS sources respectively. Further 
separation of the protostars into Class and I clumps does not yield 
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Figure 2. Distribution of clump turbulent fractions, /turb = ctnt/cs, mea- 
sured from the HARP C^^O/ = 3^2 data (see text), broken down by re- 
gion for the CLFIND (upper panels) and GAUSSCLUMPS populations (lower 
panels). 




Figure 3. Distribution of clump turbulent fractions, f[^vh = ctnt/cs, mea- 
sured from the HARP C^^O / = 3 ^ 2 data (see text), broken down by the 
different source types for the clfind (upper panels) and gaussclumps 
populations (lower panels). 



statistical differences between their average /t^rb (see Tabs. [T] and 
|2|. Unfortunately, such simple averages do not provide the whole 
picture we can see in the distributions (in Fig. [3]). For example, if 
the small but significant population of high /turb starless cores were 
eliminated it would seem that protostars have marginally higher 
/turb oil average. However, the large uncertainties on the average 
/turb do emphasize the need for larger samples of objects. For the 
CLFIND population (but not the GAUSSCLUMPS sources) most of 
the high /t^j-b protostars are Class 0, and all of the low /urb Class I, 
which is perhaps what we might expect if outflow power decreases 
from the Class to I stage (Bontemps et al. , 1996[ [Paper II| . Over- 
all, there are few significant differences between the starless and 
protostellar distributions. 

We can draw interesting comparisons with the C^^O (/ = 
1) and N2H+ /turb distributions presented by |KJT07 
|4]and|5|. Our distributions are closer to the C^^O distribu- 



2 -^ 
Figs 



tions of KJT07 than their N2H+ results. The high levels of non- 
thermal motions, shown by some of our sources, are simply not 
present in N2H+ data. A striking contrast between the C^^O dis- 
tributions is that many of our cores have high /urb (^ 3.5), which 
are not seen for either protostars or starless cores by |KJT07| In 
addition, many 'KJT07' starless cores have sub-thermal linewidths 
(/turb< !)• These differences may be partly caused by the spec- 
tr al resolu tion of our observations (Av = 0.15kms~^ compared 
^), rather than intrinsically different results. 



to 



KJT07 



s 0.05 km s~ 

Our slightly poorer resolution means we cannot probe as narrow 
linewidt hs and w e possibly blend many of the separate components, 
noted in KJTOTJs C^^O spectra, into one. Given the similarities it 
seems likely our C^^07 = 3^2 data probe the same regions as 
the / = 2 ^ 1 observations of |KJT07| n amely the outer parts of 
star-forming cores, referred to by KJT07| as their envelopes. This is 
probably because of C^^O freeze-out in the dense, cold interiors of 
star-forming cores. Such regions should see an enhancement in the 
abundance of say N2H+ relative to C^^O. KJT07 investigated this 
possibility for their similar population of candidate cores in Perseus 
by measuring the variation in the integrated C^^O-to-N2H+ ratio 
with peak SCUBA flux density (an approximate proxy for central 
H2 volume density). They found high C^^O-to-N2H+ ratios (i.e. lit- 
tle C^^O depletion) occurred mainly in starless cores, which have 
the smallest SCUBA fluxes (i.e. lowest densities), whereas high- 
flux cores (mainly protostars) have low ratios suggesting high lev- 
els of C^^O depletion. 

Given that our starless and protostellar C^^O distributions are 
quite similar, we may also conclude (like |KJT07| ) that protostars 
do not affect their environment significantly. The starless distribu- 
tion may differ from the true prestellar one as starless clumps with 
the largest linewidths are likely to be unbound, and thus the true 
prestellar turbulent fractions will be smaller. Finally, we also note 
little correspondence with the simulations of |Klessen et al.| ( [2"005| ), 
whose results for a large-scale driven (LSD) gravoturbulent model 
are a lso plotted in Figs. [4] and [5] KJT07 claim to have a good fit to 
the |Klessen et al. data with their C^^O starless population, although 
the model was designed to match N2H+ observations. 



3.2 Virial theorem 

If the C^^O data trace material that will go on to form the final star 
and not just the ambient gas in which a YSO is embedded, the core 
linewidth can inform us about its stability, i.e. whether the cores 
are gravitationally bound. We can then distinguish bound, starless 
i.e. prestellar cores from unbound ones that will dissipate without 
forming a star. 

The virial theorem is often used to examine this stability and 
we will apply it to our two populations of clumps identified in 
SCUBA 850 ^tm data. The clump virial mass, Myir, is used as an 
estimate of the internal energy of the clump whilst its dust mass at 
850 i^m, M850, is as an estimate of the potential energy. If we as- 
sume that the linewidths only reflect gravity and we have spherical 
clumps the virial theorem is 



Mvi 



aa'^R 



(3) 



where o is the average three-dimensional velocity dispersion, R the 
clump radius and a a constant that depends on the form of the den- 
sity profile (a derivation can be found in e.g. |Binney & Tremaine| 
|2008|[Rohlfs & Wilson|2004]). If we assume a power law density 



(see profile, p (r) ©c r ", then (e.g. MacLaren, Richardson & Wolfendale 



[T988l l: 
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Figure 4. Distribution of /turb. for the starless clump populations identified 
with CLFIND (left panels) or GAUSS CLUMPS (right panels), measured from 
the HARP C^^ O / = 3 -^ 2 data. Overlaid for comparison are the distribu- 
tions found by |KJT07| for 157 candidate cores using C^^O/ = 2^1 (upper 
panels, pink) and N2H+ / = 1 ^ (upper panels, blue). A prediction for 
cores from the gravoturbulent simulation of |Klessen et al.|{2005) is also 
shown for turbulence driven on large scales (lower panels, brown). For that 
simulation the turbulent fraction extends beyond 3.5 to 4.3 (not plotted), 
though there are small numbers of clumps in those bins (~ 2 per cent). 
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Figure 5. Distribution of /turb, for the protostellar clump populations iden- 
tified with CLFIND (left panels) or gaussclumps (right panels), measured 
from the HARP C^^ O/ = 3 -^ 2 data. Overlaid for comparison are the dis- 
tributions found by KJT07| for 157 candidate cores using C^^OJ = 2^ 1 
(upper panels, pink) and N2H+ / = 1 ^ (upper panels, blue). A prediction 
for cores from the gravoturbulent simulation of Klessen et al. (2005) is also 
shown for turbulence driven on large scales (lower panels, brown). For that 
simulation the turbulent fraction extends beyond 3.5 to 4.3 (not plotted), 
though there are small numbers of clumps in those bins (~ 2 per cent). 



5-2n 



-n 



(4) 



Provided the C^^O line of FWHM, Av^iso, traces the bulk of the 
gas with a Gaussian velocity distribution and the same non-thermal 
linewidth as H2, we can estimate 



81n2 



Av; 



CI80 



+ - 



81n2y^r / 1 



riQi^Q 



(5) 



where m^ is the proton mass, m = 2.33 to factor in the abundance of 
He relative to H2 and m^isQ = 30. This yields for a shallow power 
law clump, n = 1.5, consistent with the shapes found by [Enoch] 



Mvi 



'HiTj) 



Av, 



CI80 



Ikms" 



Mi^) 



R 



0.02pc 



Mr; 



(6) 



We expect objects in equipartition to have Myij. ~ M^^q whilst those 
that are self- gravitating should have Myir < 2M85o. 

We take data on the clumps' dust properties from [Curtis &| 
[Richer (2010). Each clump's 850 \mi mass (Mgso) assumes the dust 
is optically thin, has an opacity, K = 0.012 cm^ g~^ and is at a sin- 
gle temperature, Tq. Again these temperatures are taken as the NH3 
kinetic temperatures ( jRosolowsky et al.|2008) ), where available, or 
10 K (for starless cores) and 15 K (for protostars) where not. The 
core radius, /^dec is the geometric mean of the two core semi-major 
and -minor axis 'sizes' each deconvolved with the beam size. These 
'sizes' are the standard deviation of the pixel coordinates about the 
core centroid, weighted by the pixel values. 

There are considerable uncertainties in any dust and virial 
mass estimates. The errors are very difficult to quantify for individ- 
ual clumps without detailed modelling. Thus, we follow comple- 
mentary studies (e.g. [Buckle et al.|2010[ [Enoch et aL|2008j ) and do 
not attempt to account for the uncertainties in our analysis, except 
for a discussion of their magnitude, which follows. The dust masses 
depend on the assumed distance to Perseus and the dust properties 
(temperature and opacity). We try to minimize the effects of dust 
temperature by using the NH3 kinetic temperature as an estimate 
of Td, which should be an accurate measurement at high volume 
densities (> lO^cm"^, e.g. Galli, Walmsley & Gongalves 2002), 



where the gas and dust are thermally coupled. Nevertheless our dust 
temperature estimates still do not account for variations in the dust 
temperature across a clump. A range of distances have been used in 
the literature for Perseus (220 to 350 pc see e.g. Paper 1) and indeed 
it may not be a contiguous cloud at a single distance. These result 
in an uncertainty of a factor of ~ 5 in the dust mass estimates. The 
virial masses depend on the assumed clump profile, with steeper 
profiles producing smaller masses, however it only varies by a fac- 
tor of 1.7 between a constant density and 1/r^ profile. This com- 
bined with the uncertainties in the distance and linewidth produce 
a gain around a factor of ~ 5 in uncertainty. As we previously noted 
(apj]), it is likely that the C^^O / = 3 ^ 2 line only traces the 
envelope of a clump, which could lead to an over-estimate of the 
non-thermal H2 linewidth (compared to estimates from say NH3 
or N2H+) and correspondingly the virial mass. Even if our esti- 
mates of the dust and virial masses were infallible, assuming that 
a clump with Myir ^ ^850 is unbound may be misleading; exter- 
nal pressure and/or magnetic fields may contain such a clump. For 
instance |KJT07| calculate that cores with external pressures con- 
sistent with their previous Bonnor-Ebert sphere modelling (,Kirk,| 



[Johnstone & Di Francesco|2006| ), should be considered in equipar- 
tition not merely self- gravitating if Myir ^ 2M85o. 

The C^^O virial masses, which we plot in Fig. [6] against the 
850 i^m mass, lie in the range 0.8 to 22.9 M0 for the CLFIND 
sources and 0.4 to 17.8 M© for the GAUSSCLUMPS sources, re- 
flecting their smaller radii. Most of the clumps lie scattered near 
the 'equipartition' line, M^^^ = M^^O- The results are similar for 
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Figure 6. Virial mass versus clump dust mass for the CLFIND (left panels, 
(a), (c) and (e)) and GAUSSCLUMPS (right panels, (b), (d) and (f)) pop- 
ulations, differentiated by region (top panels, (a) and (b)) or the classifi- 
cations of HFR07 (the rest). Virial masses are estimated from the C^^O 
/ = 3 ^ 2 Hnewidth at the catalogue positions of Curtis & Richer (2010). 
The masses are calculated using the temperatures from Rosolowsky et aD 
( |2008) where possible or 10 and 15 K for starless and protostellar cores re- 
spectively where not. Lines denote where we expect cores to be in equipar- 
tition, i.e. Mvir = Mgso (dashed) or the limit where they are self-gravitating, 
Mvir = 2M850 (dot-dashed), although the reader should note the consider- 
able uncertainties in both masses. Only clumps which have a ,HFRQ7, source 
association are plotted. 



the different algorithms, regions and evolutionary types (see Tables 
|3]and|4]for summaries of the population statistics). Least-squares 
straight-line fitting results in poorly-constrained power law expo- 
nents of 0.8-0.9 for almost all the correlations. 

Given the uncertainties in both the virial and dust mass esti- 
mates, it is not possible to draw unequivocal conclusions about the 
stability of individiual clumps. However, we will compare distinct 
populations. The different regions show similar behaviour with a 
lot of scatter. The LI 448 population is outlying for both algorithms 
with a high M^so/^yir ratio. LI 448 is interesting as it contains 
large SCUBA flux densities and massive outflows but has relatively 
weaker C^^O emission. The NGC 1333 population has a smaller 
M85o/Mvir ratio than the others, which we might interpret as be- 
ing more unbound and therefore transient, with a larger number of 
cores above the Myir = M^^q threshold. This could be a result of its 
more perturbed environment. 

Protostellar clumps, by definition should be gravitationally 
bound. If starless clumps occupy a similar Myir — M^^q sample 
space, then this suggests that they too are a gravitationally bound 



Table 3. Summary of virial mass statistics for the CLFIND SCUBA clumps 
identified in Curtis & Richer (2010), which have a HFR07 classification. 
All the quoted errors (cr) are errors on the mean {o /\/N) not simple sample 
deviations. 



Population 


No 


Mvir 

(Mo) 


(Mo) 


MgSo/Mv/r 


erratic 


All 


50 


7.8 


0.9 


1.32 


0.13 



NGC 1333 


24 


8.9 


1.2 


0.96 


0.15 


IC348 


16 


4.0 


0.8 


1.38 


0.22 


L1448 


5 


12.6 


3.9 


2.7 


0.5 


L1455 


5 


7.9 


1.6 


1.4 


0.6 


Starless 


17 


4.6 


0.8 


1.4 


0.2 


Protostars 


33 


9.6 


1.2 


1.3 


0.2 


Class 


19 


10.1 


1.4 


1.3 


0.2 


Class I 


14 


9 


2 


1.2 


0.3 



Table 4. Summary of virial mass statistics for the GAUSSCLUMPS SCUBA 
clumps identified in jCurtis & RicherH2010} , which have a HFROTj classifi- 
cation. All the quoted errors (cr) are errors on the mean (a/y/N) not simple 
sample deviations. 



Population 


No 


Mvir 

(Mo) 


(Mo) 


MgSo/Mv/r 


erratic 


All 


64 


7.1 


0.7 


0.9 


0.1 



NGC 1333 


30 


8.9 


1.1 


IC348 


17 


5.0 


0.8 


L1448 


12 


6.5 


1.5 


L1455 


5 


4.3 


1.1 


Starless 


24 


6.0 


0.6 


Protostars 


40 


7.8 


1.0 


Class 


27 


8.4 


1.1 


Class I 


13 


7 


2 



0.68 


0.08 


0.75 


0.14 


2.2 


0.4 


0.8 


0.3 


0.87 


0.16 


0.91 


0.13 


1.06 


0.18 


0.67 


0.14 



population. In Fig. [6] the starless and protostellar clumps occupy 
the same regions of the plot and their mass ratios are very similar 
in Tabs. |3] and [4] This implies that the starless population is gravi- 
tationally bound as well and truly prestellar in nature (as found by 
[Enoch etal.pOOS i). 

The clumps also seem spread around the 'equipartition' 
line Myir = MgsQ. Approximate equipartition was also found by 



Caselli et al. (2002) in a sample of 60 starless cores using N2H^ 



Tachi- 



linewidths. Other studies e.g. the C^^O cores looked at by 

|hara et aL] ( |2002| ) found the virial mass considerably larger than 
the core mass estimate. The latter result is successfully explained 
by the gravoturbulent models of |Klessen et al.|f20 05), who find 
a large majority of starless cores with Myir ^ M for both large 
and small driving- scale turbulence models, with their protostellar 
sources non-coincident in parameter space at M > Myir. Our results 
would not favour such a model for Perseus, even in its clustered 
environment. 



3.3 Localized velocity gradients 

Systematic variations in the C^^O line centre velocity are apparent 
across the face of many of the identified SCUBA clumps, e.g. the 
velocity gradually increases along a particular direction. There are 
a number of plausible explanations for such gradients: (i) rotation, 
(ii) outflows, (iii) motions between smaller unresolved constituent 
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clumps or (iv) if gravoturbulent models are true, cores form at the 
stagnation points in convergent flows (e.g |Padoan et al.|2001| ) and 
velocity gradients may arise from colliding gas streams. It is diffi- 
cult to distinguish between each scenario so we follow the majority 
of studies and focus on analysing the rotational properties of the 
clumps from measured velocity gradients. 

A strong motivation to understand the details of rotation in 
star-forming cores is provided by the discs out of which most, if 
not all, stars are born (e.g |Shu et al.|1987] l. Planets are thought to 
originate inside such protoplanetary discs and the details of their 
formation crucially depend on various disc parameters, such as 
surface density, controlled by the detailed evolution of the angu- 
lar momentum of the parent star-forming core (e.g. |Lissauer|1993[ 
|Ruden|1999] ). Additionally, there is the classical 'angular momen- 
tum problem' of star formation (e.g. |Spitzer|1978| ): the angular mo- 
mentum of prestellar cores is orders of magnitude larger than that 
which can be contained within a single star, even though cores are 
observed to be rotating much less than originally predicted ( jGood-| 
man et al.|l993 Caselli et al.,200 2p] A plausible solution is pro- 



vided by magnetic braking at the early low-density phases, where 
the field lines are strongly coupled to the gas and transfer angular 
momentum from the contracting core to the surrounding medium 
(e.g. |Mouschovias|1987| ). Recent results suggest that gravitational 
interactions may dominate over magnetic braking. MHD models of 
self-gravitating, decaying ( [Gammie et al.|2003] ) and driven ( |Lietal.| 
■2004) turbulence have angular momenta consistent with each other 
and observations butJappsen & Klesse n (2004 hereafter [JK04| l and 
[Tilley & Pudritz|p004| find similar results in purely hydrodynamic 
frameworks. Even if magnetic fields do dominate, they cannot in- 
definitely strip away angular momentum or support a core against 
collapse because ambipolar diffusion will eventually lead to dy- 
namic collapse. Finally, for cores to fragment into multiple systems 
some angular momentum must be present and perhaps cores with 
the largest angular momenta will go on to form binaries (e.g. |Lar-| 
[son 2003 , Goodwin et al. 2007 ). 

The pre-eminent observational study in this area is the search 
for solid-body-rotation in ~40 NH3 cores by [Goodman et al.| ( |l993| 
hereafter |GBF93| ). In all their objects the rotational energy is at 
most a few per cent of the gravitational and cannot provide sup- 
port. An evolutionary sequence can be built using their results with 
others (e.g. [Caselli et al . 2002), demonstrating the specific angular 
momentum (j = J/M) decreases with decreasing scale (JK04J. 

Rotational signatures become more complicated at higher res- 
olution and once collapse has started in protostars. |Belloche et al.| 
([2002) studied a young Class protostar, IRAM 04191+1522, 
finding two distinct regimes of collapse: the inner (r < 2000 — 
4000 AU), rapidly collapsing and rotating whilst the outer (4000 < 
r < llOOOAU) has only moderate infall and rotation. The fall in 
rotational velocity beyond the 4000 AU boundary and the flat inner 
profile suggests that the inner region's angular momentum is con- 
served whilst it is dissipated in the outer perhaps by magnetic brak- 
ing. At the distance to Perseus the r ~ 4000 AU boundary is bigger 
than the JCMT beam (32arcsec diameter compared to 15arcsec) 
so we may be able to probe the inner regions, although C^^O will 



^ Furthermore, |Dib et al.| ^2010} point out that measurements of core an- 
gular momenta from global velocity gradient fitting tend to overestimate 
the intrinsic (three-dimensional) angular momenta by a factor of ~ 8- 
10, as complicated fluctuations in the three-dimensional velocity field are 
smoothed out. 



possibly freeze-out at < 5000 AU in starless and young protostellar 
cores. 



3.3.1 Velocity gradient fitting 

A clump undergoing solid-body rotation will display a linear ve- 
locity gradient across its face perpendicular to the rotation axis, we 
follow |GBF93 I and flt a linear gradient, Vvlsr. across each clump 
using: 



vlsr = vo + aAa + bAd 



(7) 



where vq is the systemic clump velocity, Aa and A5 the angular 
ascension and declination offsets from the clump centre and a and 
b the projections of the gradient per radian on to the a and d axes. 
The gradient has a magnitude, ^ = |Vvlsr| = (^^ + b^Y'^ /D, 
where D is the distance to the object, at an angle (east of north), 
% = arctan(a//7). 

We explore theC^^O/ = 3^2 velocity fleld of our two pop- 
ulations of SCUBA 850 i^m clumps. At each SCUBA map pixel 
the corresponding spectrum in the HARP data has been fitted with 
a Gaussian profile to extract its centre velocity and linewidth. Fits 
were only performed where the spectral peak was greater than three 
times the rms noise, estimated on a line-free portion of the spectrum 
for every spatial position. LI 45 5 was omitted, as its weak C^^O 
emission had too few fits. We then used the Levenberg-Marquardt 
algorithm implemented in the scientific python module, SciP3r| 
to perform non-linear least-squares fitting of Equation [T] to every 
clump about its SCUBA peak. |GBF93| performed the same analy- 
sis on sets of randomly generated maps with no gradient and report 
that approximately 10 per cent of those random maps were found to 
have significant gradients if a "3 a" criterion is used i.e. ^ ^ 3c%. 
Using the same significance criterion we expect similar levels of 
reliability. 



3.3.2 A signature of rotation ? 



We noted there are four plausible causes of velocity gradients. In 
this section, we look at some of the suggestions other than rota- 
tion in more detail. First, some of our clumps may be composed of 
smaller constituents. The gradient then might measure the velocity 
dispersion of the multiple cores (GBF93 ). If the SCUBA clumps 
are gravitationally bound, as suggested in §3.2| such dispersion in 
a multi-core clump may signal its rotation about a common centre. 
|HFR07 1 estimate within their similar core catalogue, > 10 per cent 
of clumps will break into multiple sources, by comparing the most 
luminous sources to higher resolution observations. 

Second, the exact region that is assigned to each clump can 
change the fitted gradient's magnitude, direction and level of signif- 
icance. One of |GBF93"[ s original sample, TMC-IC, a starless core 
in Taurus, has been mapped over a larger area at higher resolution 
by |Schneeet al. (2007 ). They found a more complicated pattern of 
local gradients no longer consistent with solid-body or differential 
rotation, which had previously been noted by |Caselli et al.| ( [2002| ). 

Finally, if the C^^O line centre is affected by outflows, then 
this could be misinterpreted as a sign of rotation. |GBF93T also fit- 
ted velocity gradients to C^^O 7=1^0 data in a sub-sample of 
their cores - both with and without embedded objects - and found 
they were in the same direction as the gradients from NH3. How- 
ever, they did not formally examine any outflows. Towards starless 
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Table 5. Approximate outflow orientation compared to the velocity gradient 
direction for CLFIND clumps with clear bipolar outflows. Classifications are 
from lHFR07l 



Table 6. Number of clumps with significant velocity gradients by region, 
source type and algorithm. Percentages give the detection rate for the par- 
ticular type of source in each region or overall. 



Sub-region 



CLFIND 

ID 



Class 



(degEofN) 



Outflow orientation 
(degEofN) 



Source Type 



NGC 1333 


6 





IC348 


1 





IC348 


2 





IC348 


4 





L1448 


1 





L1448 


2 





L1448 


4 






126 

-51 

50 

4 

139 
78 
154 



-113 

-64 

68 



124 
165 
168 



cores, outflows should not be a problem. To explore the effect on 
our gradients, we compare the direction of the velocity gradient to 
the orientation of the outflow from every protostellar source with 
a significant gradient in the CLFIND catalogue. We estimated the 
outflow position angle from our CO / = 3 ^ 2 datacubes (a full 
outflow analysis was presented in [Paper Il] l. For some of the proto- 
stellar clumps, the outflows are complicated and we have ignored 
such cases, only comparing clumps where a characteristic bipolar 
structure is apparent (see Table [5]). In a number of these cases the 
velocity gradient is in the same direction as the outflow, implying 
we are tracing the outflow rather than solid-body rotation with the 
gradient. If we were viewing solid-body rotation of the core or a 
disc around it, we would expect the gradient and outflow to lie per- 
pendicular to one another. 

The data themselves present a more complicated picture. In 
Fig.]?] we overlay the outflows on top of the C^^O centre veloci- 
ties for the sources in Table pi Often the higher C^^O line centres 
(i.e. red-shifted gas) correspond closely to the extent of the red- 
shifted CO outflow lobe and similarly for the blue- shifted gas. For 
instance, the complex outflow structure in L1448-1 is closely fol- 
lowed by the line centre data. Furthermore, the blue-shifted outflow 
lobe in NGC 1333-6 is coincident with blue-shifted C^^O line cen- 
tres. The outflow itself extends away from the clump to the west, 
resulting in orientations that seem very disparate in Table |5] In 
two cores there seems to be little outflow-gradient correspondence: 
IC348-2 and L1448-4. Perhaps it is more convincing to argue that 
none of the cores exhibit gradients perpendicular to their outflows. 
Therefore, there appears to be a link between the outflow orienta- 
tion and the velocity gradient, although how much this influences 
the gradient found in every protostar is difficult to quantify. 

Thus, we must view the interpretation of a linear velocity gra- 
dient as unequivocal evidence of solid-body rotation with some 
scepticism. However, in the absence of a strong outflow i.e. for star- 
less clumps and assuming we can consider each clump as a single 
object, rotation would seem to be a likely cause of a coherent veloc- 
ity gradient and the only one we consider in the following analysis. 



3.3.3 Fitted gradients 

Table [6] lists the rates of detecting significant velocity gradients to- 
wards various core populations, separated by their different clump- 
finding algorithms, regions or types of source. They are similar 
for both algorithms with approximately a third of clumps display- 
ing significant velocity gradients. The detection rate for Class 
sources in the CLFIND catalogue is nearly twice as large as for the 
other types and Class Os with GAUSS CLUMPS. This is probably be- 
cause Class sources have typically stronger outflows than Class Is 



NGC 1333 



Region 
IC348 



L1448 



Total 



CLFIND Population 



Starless 
Class 
Class I 

No 



HFR07 



All 



ID 



4(67 %) 

7(64 %) 

4(44 %) 

1(8 %) 

16(41 %) 



2(17 %) 
3(100 %) 
1(50 %) 
1(17 %) 
7(30 %) 



0(0 %) 

3(75 %) 

0(0 %) 

1(9 %) 

4(25 %) 



6(32 %) 

13(72 %) 

5(45 %) 

3(11 %) 

27(35 %) 



GAUSSCLUMPS Population 
Starless 4(50%) 6(55%) 0(0%) 

Class 6(50%) 1(33%) 0(0%) 

Class I 3(33%) 0(0%) 0(0%) 

No HFR07I D 11(31 %) 3(38 %) 1(6 %) 
All' ' 24(37%) 10(42%) 1(4%) 



10(48 %) 

7(35 %) 

3(35 %) 

15(24 %) 

35(31 %) 



(e.g. |Bontemps et al.|1 996, Paper 11). It is clear in Fig.|7]that many 
of the outflows extend beyond their associated clump. However, 
the clumps found by GAUSSCLUMPS are smaller ( [Curtis & Richerj 
|2010| ) with a definite elliptical shape, which does not necessarily 
follow the orientation of the outflow. The lower detection rate for 
Class sources with GAUSSCLUMPS is perhaps another example 
of how their velocity gradients are dominated by outflows. If the 
smaller GAUSSCLUMPS cores trace their shape less effectively, we 
would expect a lower detection rate. 

The fitted gradients are in the range 1.5 to 13.0km s~^ pc~^ 
and 2.5 to 16.0kms~^ pc~^ for CLFIND and GAUSSCLUMPS re- 
spectively. We provide a summary of all the different rotational 
properties for different categories of clumps in Tables |7]and|8] Fig. 
|8]is a histogram of their magnitude by region. The distributions are 
qualitatively similar for the different regions with the usual uncer- 
tainty due to the small sample sizes; Kolmogorov-Smirnov (K-S) 
tests do not conclusively reject or confirm the hypothesis that each 
sample is drawn from the same population. On average across all 
the regions, the gradients are (^) = (5.7 ± 0.5) kms~^ pc~^ for 
CLFIND and (^) = (6.9±0.6)kms~^pc~^ for GAUSSCLUMPS. 
The fitted GAUSSCLUMPS gradients are larger than those for 
CLFIND but may be misleading as there are only 16 clumps with 
detections in both. We plot the gradient found with the CLFIND des- 
ignation versus that with GAUSSCLUMPS for these common sources 
in Fig. |9] There is fairly good agreement on the gradient with the 
protostellar sources deviating the most from the line of equal gradi- 
ents. This might again be because the C^^O line centres are affected 
by outflowing gas which is not well traced by the GAUSSCLUMPS 
outlines. 

There are too few sources of each classification to draw firm 
conclusions about any trends with age. It does seem in the distribu- 

lOkms-^ 



tions of Fig. 



10 



that the starless cores have 



pc 



and larger gradients are protostellar. The protostellar average is 
larger than the starless one for the GAUSSCLUMPS population: 
(^) = (7.0 ±1.0) compared to (^) = (5.0 ±0.6) km s-^ pc-^ Al- 
though they are similar for CLFIND (^) = (5.8 ±0.7) compared 
to (^) = (5.5 ± 1.0) kms~^ pc~^ This might point to an outflow 
contribution to the gradient. 

The gradients we find i n Perseus are much gre ater than those 



GBF93 



Casellietal.| ( |2002| ) in their N2H+ 



originally found by 

/ = 1 ^ survey of 60 low mass cores: 0.3 to 3.9 and 0.5 to 

6.0kms~^pc~^ respectively. Outflows cannot explain all of this 
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Figure 7. Outflow orientation compared to the velocity gradient direction for clumps in Tablejs] The colour-scale is the C^^O / = 3 ^ 2 line centre in km s ^ 
for map pixels allocated to each clump. Overlaid are contours of CO / = 3 ^ integrated antenna temperature / T^dv for the blue- and red-shifted line wings, 
integrated from: -5 to 3 (blue) and 12 to ISkms"^ (red) for NGC 1333, -5 to 6 (blue) and 11 to ISkms"^ (red) for IC348 and -5 to (blue) and 8 to 
15kms~^ (red) for L1448. Contours are at various heights in Kkms~^ (n is an integer from zero to ^max): (a) (2-1-4^), ^max = 4 (blue) and (4+4^), /imax = 6 
(red), (b) (6+2n), ^max = 6 (blue) and (6+2n), rirnax = 3 (red), (c) (^+l.5n), ^max = 6 (blue) and (4+l.5n), n^ax = 8 (red), (d) (5+n), n^ax = 4 (blue) and (2+n), 
^max = 2 (red), (e) i2+n), n^ax = 8 (blue) and (2-1-2^), ^max = V (red), (f) (2+4/2), ^max = 5 (both blue and red), (g) (2-i-4«), n^ax = 5 (blue and red). 



Table 7. Average rotational properties derived from the velocity-gradient fitting for the CLFIND clump population. All the errors listed (cr) are errors on the 
mean and not sample deviations. 



Population 


Number 


^ 


(7^ 


ftot 


^/3 


J 


^j 






(kms~ 


^pc-1) 






(kms" 


"^pc) 


All 


27 


5.7 


0.5 


0.014 


0.003 


0.0019 


0.0002 


NGC 1333 


16 


5.7 


0.7 


0.0040 


0.0010 


0.0017 


0.0002 


IC348 


7 


5.0 


0.6 


0.015 


0.006 


0.0017 


0.0003 


L1448 


4 


7 


2 


0.032 


0.016 


0.0032 


0.0009 


Starless 


6 


5.5 


1.0 


0.013 


0.007 


0.0021 


0.0005 


Protostars 


19 


5.8 


0.7 


0.008 


0.003 


0.0019 


0.0003 


Class 


14 


6.1 


0.8 


0.010 


0.004 


0.0021 


0.0003 


Class I 


5 


4.9 


1.5 


0.003 


0.001 


0.0011 


0.0001 
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Table 8. Average rotational properties derived from the velocity-gradient fitting for the GAUSSCLUMPS clump population. All the errors listed (cr) are errors 
on the mean and not sample deviations. 



Population 


Number 


(kms~ 


^pc-i) 


ftot 


^/3 


i 

(kms" 


-ipc) 


All 

NGC1333 

IC348 


35 
24 
10 


6.9 

7.4 
5.6 


0.6 
0.7 
0.8 


0.013 
0.015 
0.007 


0.004 
0.006 
0.001 


0.00120 

0.0012 

0.0011 


0.00014 

0.0002 

0.0002 



Starless 


10 


5.0 


0.6 


0.006 


0.001 


0.0010 


0.0002 


Protostars 


10 


7.0 


1.0 


0.009 


0.004 


0.0012 


0.0002 


Class 


7 


6.4 


1.0 


0.008 


0.005 


0.0012 


0.0002 


Class I 


3 


8 


3 


0.012 


0.010 


0.0012 


0.0003 
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Figure 8. Distribution of significant (cr^/^ ^ 3) velocity gradients for the 
CLFIND (top) and GAUSSCLUMPS (bottom) clump populations. The L1448 
distribution is not plotted for GAUSSCLUMPS as it has only one significant 
gradient. 




Figure 10. Distribution of significant (cr^/^ ^ 3) velocity gradients for the 
classifications of ,HFR07, forthe CLFIND (top) and GAUSSCLUMPS (bottom) 
clump populations. 




Figure 9. Fitted velocity gradient for GAUSSCLUMPS versus CLFIND des- 
ignations for the 16 clumps which have significant detections with both 
algorithms, by region (left) or by the classifications of |HFR07| (right). The 
line marks where both gradients are equal. 



difference because our starless clumps have significantly larger gra- 
dients as well. The discrepancy is reminiscent of that found by 
GBF93|towards B 361, which was also investiga ted by | Arquilla &| 



Goldsmith| ( [l985| with ^^CO 7=1^0 data. |GBF93| estimated 



a gradient five times smaller t han the earlier study, 0.7 ver sus 



3.3kms~^ pc~\ which had led Arquilla & Goldsmith ( 1985i to 



|GBF93 'explain the difference by emphasizing that tracers of lower 
critical densities can produce complicated velocity patterns dom- 
inated by outflows or clump-to-clump motions that mimic solid 
body rotation. However, |01mi, Testi & Sargent] ( [2005| ) fitted gra- 
dients to a number of cores in Perseus, finding similar magnitudes 
for the two common cores also in |GBF93] and generally in the same 
range 0.14 to 2.32 kms~^ pc"^ Their study used C^^O 7=1^0, 
CS 7 = 2 ^ 1 and N2H+ 7=1^0 data with their C^^O fit- 
ting not displaying markedly different gradients. A critical factor 
will undoubtedly be resolution. The angular resolution of previ- 
ous surveys is worse than the 17.7 arcsec of our maps: ~88, ~54 
and -46 arcsec for GBF931 [ Caselli et al.| ( [20 02]) and |01mi et aT] 
(2005 ) respectively. Hence even for their closest sources in Taurus 
(at 140 pc), their best linear resolution (0.037 pc) is nearly twice 
as large as ours (0.021 pc). A larger beam will tend to smooth out 
differences between regions and reduce the overall gradient mag- 
nitude. This effect may be large enough to explain the differences 
for cores perhaps three times as distant as Perseu^ in a beam four 
times as large as ours. 



conclude that rotation was dynamically significant in dark clouds. 



^ The most distant sources examined by jCaselli et al, 
at 900 pc and L1389 at 600 pc. 
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Figure 11. Velocity gradient versus axis ratio for CLFIND (top panels) and 
GAUSSCLUMPS sources (bottom panels) broken down by region (left pan- 
els) or by the classifications of |HFR07| (right panels). 
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Figure 12. Distribution of the angle between the velocity gradient and the 
major axis of the GAUSSCLUMPS clumps by region (left) or by the classi- 
fications of HFR07 (right). If rotational motions are elongating the clumps 
we would expect the gradient and axis to be in the same direction making 
the angle zero. 



3.3.4 Gradient orientation 

Centrifugal stresses due to energetic rotation will flatten cores 
along the rotation axis. This tends to produced oblate cores 
(|GB F93'), currently the favoured shape (e.g. Jones, Basu & Dubin-| 
|ski| 2001), which could also be caused by strong magnetic fields. 
If rotation dominates, we might expect the the degree of core elon- 
gation to increase with increasing velocity gradient and a core's 
major axis to lie parallel to the gradient i.e. perpendicular to the 
rotation axis. Fig. pT] shows no relationship between the magnitude 
of the velocity gradient and the elongation of a clump, quantified 
through the axis ratio (taken from jCurtis & Richer|2010| ). It comes 
as no surprise then that the angle between the clump major axis 
and the velocity gradient in Fig.[T2](only computed for the GAUSS- 
CLUMPS sources), seems to be distributed at random as well. There- 
fore, clump rotation is unlikely to be energetically significant. 

Is there a correlation between the orientation of the rotation 
axes of neighbouring clumps or an anti-correlation to keep the net 
angular momentum small? In Fig. [13] we overlay the velocity gra- 
dients on dust emission maps of the various regions. Reassuringly 
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Figure 13. Spatial distribution of the velocity gradient detections in 
NGC 1333 (top), IC348/HH211 (middle) and L1448 (bottom). Arrows are 
drawn from the peak position of the associated clump for clumps where 
there are significant detections in either the CLFIND or GAUSSCLUMPS cat- 
alogues. The direction and length of the arrow represents the direction and 
magnitude of the gradient and the detections are labelled with the reference 
from our SCUBA catalogues (Curtis & Richer 2010). The colour-scale is 
SCUBA 850 |Lim emission from to 1600 mJy beam"^ . 
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the directions and magnitudes of the gradients measured using both 
algorithms' allocations match very well. In small patches the ar- 
rows of neighbouring clumps line up but on large scales there seems 
to be little correlation. This is similar to the gravoturbulent mod- 
els of |JK04| where shortly after collapse starts, the angular mo- 
menta in small regions are well aligned over a moderate correlation 
length with cores further away spinning in random directions. Their 
explanation is that neighbouring cores are formed from the same 
reservoir of material and therefore can be expected to have similar 
angular momenta. During subsequent accretion and evolution, the 
associated correlation length decreases and neighbours lose their 
alignment, as embedded cores are ejected and turbulence disrupts 
the current material or brings in new gas. 




Figure 14. Distribution of the ratio of rotational to gravitational energies, 
jSrot- For comparison the distributions found by .GBF93 (green) and (Caselli] 
|et al.|2002| orange) are overlaid. 



33.5 Dynamic support? 

To quantify the level of support that the velocity gradients might 
provide against gravitational collapse we calculate jSrot, the ratio 
of rotational to gravitational energy. [GBF93 1 define this ratio for a 
uniform density sphere: 



fto 



p3^2 



(1/2)70)^ 

3GM^/5R~ SGMsin^/ 



(8) 



where the moment of inertia is / = pMR^ with p = 2/5 for a uni- 
form density sphere and the angular momentum is co = ^ / sin i with 
/ the angle of inclination to the line of sight. Assuming values rep- 
resentative of these clumps and sin/ = 1, this becomes: 



ftot = 6.2 X 10" 



0.02 pcy \Mq 



1 km s 1 pc 



(9) 



Across all the clumps the average values of jSrot are: (/3rot) = 
(0.014 ± 0.003) and (0.013 ± 0.004) for CLFIND and GAUSS- 
CLUMPS respectively. At most the rotational energy is just six per 
cent for the CLFIND sample and 14 per cent for GAUSSCLUMPS of 
the gravitational energy in the clumps. Fig. [14] shows the distribu- 
tion of jSrot for the two clump populations. There is a strong pref- 
erence for low ratios, with the fraction of clumps rapidly falling 
off as jSrot increases. Additionally, both populations closely match 
each other and the ratios found by |Caselli et al.| ( [2'002| |, while the 
|GBF93 1 distribution is wider, spreading to higher values. K-S tests 
could not confirm the samples are drawn from the same population 
as each other or either of the examples from the literature to any 
degree of significance. Even though we found larger velocity gra- 
dients than the previous studies, jSrot is not any higher. This is prob- 
ably because our clumps are significantly smaller: on average the 
clump radii are 0.03 and 0.02 pc for CLFIND and GAUSSCLUMPS 
respectively while the , Caselli et al.| ( [2002| ) cores have an average 
of 0.06 pc. As jSrot depends on R^ this will reduce our values by 
^9 relative to theirs for an equally massive source, compensating 
for our larger derived gradients. Thus, rotation is not dynamically 
significant in star-forming cores and should not support the core 
against collapse. 

It would be interesting to examine how the level of support 
varies region-to-region and as a clump evolves. However, the num- 
bers here are too small to draw any firm conclusions, with one 
or two objects heavily distorting overall averages. This is further 
compounded by the outflow ambiguities of the protostellar sources. 
Nevertheless, there are no vast differences in the populations by re- 
gion or age with almost all the clumps having /3rot ^ 0.02. 
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Figure 15. Distribution of j by region and algorithm: CLFIND on the left 
and GAUSSCLUMPS on the right. 



3.3.6 Angular momenta 

The core rotation speed can be quantified via the specific angular 
momentum, j = J/M, i.e. the angular momentum per unit mass: 

2 

where co = ^/ sin / is again the angular velocity of the clump, R 
its radius and the factor 2/5 is for a constant density sphere. This 
yields for representative values: 

2 



:COR'^ 



(10) 



7 = 1.6x10-4 



1 km s ^ pc 



0.02 pcy 



kms pc 



(11) 



The distribution of j has been important, with modellers using the 
one derived by |GBF93 [either as an input to their models or a target 
(e.g. |Burkert & Bodenheimer|2000l|JK04| ). Furthermore, the distri- 
bution of angular momentum across a core and its evolution may 
prove whether magnetic fields are necessary to strip off angular 
momentum as collapse ensues. 

In Fig. [15] we plot the distribution of j by region. There 
is little difference between the different algorithms or clumps 
from one region to the next. For CLFIND the overall average is 
{j) = (1.9±0.2) X 10-3 kms-i pc whilst {j) = (1.7±0.2) x 10-^ 
(1.7±0.3) X 10-3 and (3.2±0.9) x 10-^ inNGC 1333, IC348 and 
L1448 respectively. The GAUSSCLUMPS data are skewed to lower 
j due to their smaller clump radii but there are still few differences 
across the regions: (j) = (1.20 ±0.14) x 10-^, (1.2 ±0.2) x 10-^, 
(1.1 ±0.2) X 10-3 kms-i pc overall and in NGC 1333 and IC348 
respectively. 

I JK04| distinguish between starless and protostellar sources in 
their gravoturbulent simulations which focus particularly on rota- 
tional properties. They find wide continuous distributions for y, 
stretching over two orders of magnitude, as we do. Their prestel- 
lar cores have an order of magnitude larger j than the protostars. 
Fig. [16] compares our distributions of j (overall and for the starless 
clumps) to previously observed distributions and |JK04| To match a 
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Figure 16. Distribution of specific angul ar momenta, j for starless cor es, (a) 
N2H+ / = 1 ^ cores of |Caseni et'ar| ( |2002, grey) compared to tlie |JK04| 
model (brown) scaled by ^ = 10^ cm~^ and T = 10 K. (b) All our clfind 
detections (red) and just those in starless clumps (black), (c) NH3 cores 
of GBF93| [Blanco & Goodinanl fT998) and | Jijina et "H!] fT^99) (green) 
compared to the |JK04| model (brown) scaled by w = 10^ cm~^ and T = 
10 K. (d) all our GAUSSCLUMPS detections (blue) and just t hose in starless 
clumps (black), (e) N2H+ / = 1 ^ high-mass cores of |Pirogov et al.| 
( 2003 > compared to the JK04 model (brown) scaled by n = 10^ c ~^ and 
r = 50K. 



particular set of observations, pKO?] scale their code output accord- 
ing to the mean density, n, and temperature, T, with j depending 
on these as j ©c T j ^. With appropriate scalings their data match 
three different observational r egimes: low-m ass cores mapped in 
the high-density tracer N2H+ ( Caselli et alJ 2002), low-mass cores 
in the slightly lower density tracer NH3 (GBF93 with Barr anco &| 
|Goodman|1998| and [Jijina, Myers & Adams, 1999j and high -mass 
cores in the high-density tracer N2H+ ( Pirogov et al.|2003] ). Our 
observed distribution has no overlap with the high-mass cores of 
[Pirogov et al.| ([2003) and little with [GBF93[ having much lower 
]. For the starless clumps only, we find (j) = (2.1 ±0.5) x 10~^ 
and (1.0 ±0.2) x 10~^ km s~ ^ pc for CLFIND and GAUSSCLUMPS 
respectively. The most similar distribution is that of [Caselli et al.[ 
( |2002| ), which is over twice as wide with a much larger average, 
(j) = 2.2 X 10~^ km s~^ pc, but does overlap with ours at low 7. For 
our average to match that of [Caselli et al.[ ( [2002| , presumably our 
clumps would have to be co nsiderably denser and/or colder tha^ 
the values assumed by JK04 oin = 10^ cm~^ and T = 10 K. In 
we noted that the critical density is expected to be /icrlt = 10^ cm 



However, in radiative transfer models of a free-falling protostellar 
envelope, we found peak H2 number densities of greater than a few 
times 10^ cm~^ were required to reproduce the C^^O/ = 3^2 
line strengths seen in our data ( Curtis 2009 ). Alternatively or addi- 
tionally, we are seeing considerably smaller angular momenta then 
those measured by previous authors. The narrow width of the dis- 
tribution compared to previous work is likely to reflect the uniform 
environment of our cores in a single cloud rather then spread over 
many different ones as in Caselli et aL[ ( [2002] ) and |GBF93, 

On the other hand, the smaller angular momenta observed in 
our starless clumps may not indicate lower degrees of rotation in 
the clump core but may reflect a more slowly rotating outer enve- 
lope, making c omparisons with tracers such as NH3 inappropriate. 
(2004) examined the asymmetric HCO^ line pro- 



Redman et al. 



files across L1689B, a prestellar core in Ophiuchus, which they 
modelled as a rapidly -rotating inner core contained within a static 
envelope. Again, the key diagnostics are the densities and temper- 
atures traced by the C^^O / = 3 ^ 2 line, which can only be prop- 
erly answered with radiative transfer modelling. We would expect 
to probe high densities, /icrit ^ lO^cm"^, but earlier we saw the 
similar non-thermal linewidths for the / = 3 ^ 2 line compared to 
the / = 2 ^ 1, suggesting we might not be detecting very different 
material. A further consideration in the low temperature environ- 
ments of starless cores is CO depletion. Estimating the degree of 
depletion in these cores is difficult without similar observations 
in "late-depleting" molecules, such as N2H+. However, models 
(e.g. Walmsley, Flower & Pineau des For ets|2004] ) suggest that CO 
freezes-out in the central 7000 AU of a prestellar core (this is really 
a tem perature and density dependence). Furthermore, we noted in 
J3.l| that there is some existing evidence for considerable C^^O de- 
pletion in candidate star-forming cores in Perseus (KJT07 ) Thus, 
our C^^O signal in evolved starless clumps and young protostars 
at least should have little contribution from the dense inner core 
and this reasoning supports the idea we are probing a more slowly 
rotating envelope. 

In JK04 protostars have a distribution of j that narrows with 
time, staying around the same average. Their distribution is also 
similar to that observed for binaries (see Fig. [Tv]). Our distribu- 
tions for all the sources are as narrow as those for the proto- 
stars alone in lJK04[ s models. For CLFIND there is a trend of de- 
creasing j with protostellar age: (j) = (2.1 ±0.3) x 10~^ and 
(1.1 ±0.1) X 10~^kms~^ pc for Class and I clumps, which is not 
reproduced with GAUSSCLUMPS: (j) = (1.2 ±0.2) x 10~^ (Class 
0) and (1.2 ±0.3) x 10"^kms"^pc (Class I). Of course in these 
sources there is potential confusion of the rotational signatures with 
outflows which is perhaps why the expected decrease in j from star- 
less to protostellar clumps is not observed. 



3.3.7 Scalings with clump size 



[GBF93[ found that many of the derived rotation parameters scale 
with radius, including ^ (oc R-^-"^) and j (oc R^-^) but not jSrot- We 
plot the same parameters as a function of clump size in Fig. [18] 
Most of the relations found by GBF93 hold for these clumps, al- 
though there is considerable scatter. GBF93 show that the relation 
for the gradient is implied from a line width- size relation (e.g. [Lar-| 
son 1981 ), Av = R^-^, for a core in virial equilibrium, we find j3rot is 
independent of R as well, implying M oc R^ for solid body rotation 
(CO = constant) or for differential ro tation (co qc R~^) that M oc R, 
The SCUBA clumps follow M oc R^ ( jCurtis & Riche r '20 10*), so we 
can infer a more complicated form of rotation than pure solid body 
or differential to ensure jSrot is independent of size. 
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Figure 17. Distribution of the specific angular momenta, j, for protostars. (a) to (d) show|JK04|s model M6k2 with different star formation efficiencies: (a) 
15, (b) 30, (c) 45 and (d) 60 per cent, (e) Binaries among nearby G stars of Duquennoy & Mayor (1991 1, (f) Binaries in Taurus (Simon 1992). (g) and (h) show 
the distributions from paper for the CLFIND and GAUSSCLUMPS designations respectively, decomposed according to the designations of HFR07, 
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Figure 18. Rotational parameters as functions of clump size for clfind 
(left) and GAUSSCLUMPS (right). ^ versus clump radius; ftot versus /^dec 
(middle panels); j versus /?dec (bottom panels). Lines mark the |GBF93] re- 
lations: ^ oc R-^-^ and j oc R^-^ . 



The scaling of j with size is of particular interest. [Ohashi et aT] 

( |l997| ) looked at rotation in protostars in Taurus using C^^O / = 
1^0 interferometric spectra. They found at distances < 3500 AU, 
that j was at a constant value of ~10~^ km s~^ pc. These findings 
are consistent with magnetically-controlled collapse elucidated by 
[Belloche et al.| { [2002| ), where in a dense inner core, the angular 
momentum is locked at a minimum value. This region could be in- 
terpreted as a magnetically supercritical core, decoupling from its 



subcritical environment, where j scales with R. In our data a hori- 
zontal line is almost as convincing as the slope of |GBF93] indeed 
the spread of clump sizes is not really wide enough to constrain 
the trend unambiguously. Our results are therefore consistent with 
this magnetically-controlled scenario: there is a lot of spread about 
j ~10~^ km s~^ pc but no j much below this value. How the rota- 
tional velocity changes across a clump is less clear and is a key dis- 
criminator between different theories. Such measurements would 
require high-resolution, interferometric observations of many cores 
in different stages of collapse. 



4 SUMMARY 

We have analysed the kinematic signatures of two populations of 
dust continuum clumps identified in SCUBA 850 jam data using 
either CLFIND or GAUSSCLUMPS (C urtis & Riche r 2010), across 
four regions (NGC 1333, IC348/HH211, L1448 and L1455) of ac- 
tive star formation in the Perseus molecular cloud. The kinematics 
were de rived fro m our large-scale (600arcmin^) C^^O 7 = 3^2 
survey |Paper l| ). First, the non-thermal contribution to the C^^O 
linewidth was calculated and compared to the results of |KJT07| and 
the gravoturbulent models of jKlessen et al.| ( [2005| ): 

(i) Most clumps have supersonic non-thermal linewidths, 
(cjnt/cs) = 1.76 ±0.09 (CLFIND population) or 1.71 ± 0.05 
(GAUSSCLUMPS), with overall distributions similar to the C^^O 
/ = 2 ^ 1 data of |KJT07 [ towards dense core candidates in Perseus. 
This implies we are also probing the 'envelopes' of the star-forming 
clumps rather than deep within their interiors. 

(ii) There is little difference in the level of the non-thermal 
linewidth contribution between protostars and starless clumps, im- 
plying that protostars do no affect their environment significantly. 
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Clumps in NGC 1333 and IC348 are in the vicinity of young stellar 
clusters, explaining the wider linewidths in NGC 1333 ((ont/cs) = 
1.91 ±0.07) than L1448 (1.52±0.13) for the GAUSSCLUMPS ob- 
jects (although their linewidths are similar for CLFIND population: 
2.05 ± 0. 12 and 1 .8 ± 0.2 respectively). IC348 however, has the nar- 
rowest lines of all (1.36 ± 0.09 with CLFIND and with 1.44 ± 0.07 
GAUSSCLUMPS), possibly related to its statistically younger popu- 
lation ( [Hatchell et al.|2005] ), containing more starless clumps. 

Second, the C^^O linewidths allowed an estimate of each 
clump's virial mass: 

(iii) Most clumps in Perseus appear bound and close to equipar- 
tition. The ratio of the SCUBA to virial mass for clumps based on 
flux allocation using CLFIND is (Mgso/Mvir) = 1.32 ±0.13, com- 
pared to 0.9 ±0.1 for the GAUSSCLUMPS decomposition. 

(iv) Starless clumps occupy a similar section of MgsQ-Mvir pa- 
rameter space to the protostars suggesting that they are gravitation- 
ally bound and therefore truly prestellar. This has implications for 
our previous work investigating the starless CMF in Perseus ( [Curtis | 
|&Richer2010). 



These results are in contrast to the models of IKlessen et al.l 
( |2005| ) where a majority of starless core have Myir ^ M and occupy 
a different parameter space to the protostars. 

Finally, we fitted a linear velocity gradient to the C^^O line 
centres, across the face of each SCUBA clump, again identified 
either with CLFIND or GAUSSCLUMPS. Significant detections were 
found in approximately a third of the clumps, which we interpreted 
in terms of solid-body rotation. 

(v) A correlation between CO J = 3 ^2 outflows and the C^^O 
line centres is observed, implying any gradients in outflow sources, 
i.e. protostars, may reflect outflows rather than rotation. 

(vi) The fitted gradients, in the range ~ 1 to 16 km s~^ pc~^ are 
larger than those previously observed. This probably results from 
the higher angular resolution of our observations and/or outflow 
contamination. 

(vii) There is no correlation between the gradient direction and 
the clump orientation nor between the gradient and clump axis ra- 
tio, implying that the rotation is not energetically significant. Fur- 
thermore, the ratio of the clump rotational to gravitational energy 
is typically < 0.02, demonstrating that rotation is not dynamically 
important as well, with a d istribution very similar to that of the 
N2H+ cores of |Caselli et"ari|2002 ). 

(viii) The distribution of specific angular momentum is nar- 
rower and centred around lower values, j ^ 10~^ kms~^ pc, than 
previous studies. Interpreting the results of |JK04| this suggests a 
dense r and/or colder env ironment for the Perseus clumps than seen 

10^ cm- 



inthe 
T 



Caselli et al. 



([2002^ sample (taken to have n ' 



and 



10 K). This would seem unlikely given the large linewidths we 
observed towards the clumps, although radiative transfer modelling 
does suggest we are probing moderately high densities - peak H2 
number densities of a few x 10^ cm~^ for a free-falling protostellar 
envelope ( [Curtis |20"09| . Thus, we are probably seeing lower levels 
of rotation in our clumps. 

(ix) There are no strong trends in the rotational parameters with 
radius. 



The somewhat inconclusive findings of this paper motivate 
further work on the star-forming cores in Perseus. First, the con- 
struction of accurate radiative transfer models to determine pre- 
cisely the conditions traced by the CO / = 3 ^ 2 transitions. 
Second, higher-resolution or higher-tracer-density observations. 7 



found in our clumps, suggestively hovers around t he bac kground 
value, 10~^ kms~^ pc~^ found by Ohashi et al. ( 1997| ) at radii 
< 0.03 pc in Taurus's protostellar cores. To see how exactly j varies 
with radius we need to probe across the clumps themselves, using 
higher-density tracers or higher-resolution observations. 
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APPENDIX A: DERIVED CORE DYNAMICAL AND 
ROTATIONAL PROPERTIES 

We list linewidths, virial masses and rotational parameters along- 
side ot her p ropertie s for eve ry cl ump examined in Tables |A1| |A2| 



A3 



and 



A4 



In Figs. 



Al 



and 



A2 



we plot the C^^O / = 3 ^ 2 line 



centre velocity and best-fitting velocity gradient for every SCUBA 
clump with a significant detection. 
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Table Al. Dynamical properties of SCUBA clumps found with CLFIND. The full version of this table is available as Supporting Information to the online 
version of this article. 



Sub-region ^ 


Clump ^ 


Hatchell ' 


7kH3^ 


Vc(Ci8o)^ 


Avciso ^ 


ctnt^ 


Mvir' 




ID 


ID 


(K) 


(kms-i) 


(kms-i) 


(kms"^) 


(Mo) 


NGC 1333 


1 


42 


13.5 


7.36 


1.91 


0.81 


11.3 


NGC 1333 


2 


43 


16.3 


8.24 


1.53 


0.64 


11.9 


NGC 1333 


3 


44 


16.5 


7.62 


1.19 


0.50 


12.3 


NGC 1333 


4 


41 


15.0 


7.48 


1.81 


0.77 


13.7 


NGC 1333 


5 


45 


16.4 


7.79 


1.40 


0.59 


22.9 



^ Name of sub-region map. 

^ Clump number from Curtis & Richer (2010). 

' H07 identifier. 

^ Rosolowsky et al. (2008) NH3 kinetic temperature where they exist or 10.0, 15.0 and 12.0 K for starless clumps, protostars and clumps with no H07 

identification where not respectively. 

^'^ Measured C^^O / = 3 ^ 2 line centre velocity, Vc and FWHM, Av^-^isq from a Gaussian fit to the line profile at the clump peak. 

^ Non-thermal contribution to the linewidth estimated using the temperatures of column (d). 

^^ Virial masses calculated using the temperatures of column (d). 

Table A2. Dynamical properties of SCUBA clumps found with GAUSS CLUMPS. The full version of this table is available as Supporting Information to the 
online version of this article. 



Sub-region ^ 



Clump ^ 
ID 



Hatchell ' 
ID 



(K) 



Vc(Ci8o)^ 
(kms-i) 



Av(-i8o 
(kms-i) 



(kms~ 



(Mo) 



NGC 1333 


1 


41,42 


15.0 


6.77 


1.04 


0.44 


11.9 


NGC 1333 


2 


43 


16.3 


8.18 


1.51 


0.64 


5.8 


NGC 1333 


3 


- 


16.5 


7.62 


1.19 


0.50 




NGC 1333 


4 


45 


16.4 


7.82 


1.35 


0.57 


8.7 


NGC 1333 


5 


46 


14.3 


8.40 


1.27 


0.54 


4.0 



^ Name of sub-region map. 

^ Clump number from Curtis & Richer (2010). 

' H07 identifier. 

^ Rosolowsky et al. (2008) NH3 kinetic temperature where they exist or 10.0, 15.0 and 12.0 K for starless clumps, protostars and clumps with no H07 

identification where not respectively. 

^'^ Measured C^^O / = 3 ^ 2 line centre velocity, Vc and FWHM, Av^-^isq from a Gaussian fit to the line profile at the clump peak. 

^ Non-thermal contribution to the linewidth estimated using the temperatures of column (d). 

^ Virial masses calculated using the temperatures of column (d). 
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Figure Al. SCUBA 850 \ym clumps identified with CLFIND containing significant velocity gradients, cr^/^ ^ 3 in NGC 1333, IC348 and L1448. The colour- 
scale denotes theC^^O/ = 3^2 line centre velocity, linearly scaled from the clump minimum (blue/black) to maximum (red/white). Contours of 850 /im 
flux density are overlaid (black) at 100, 200, 400, 800, 1600 and 3200 mJy beam" ^. A white/black circle marks the position of the clump peak emission with 
a white line positioned to show the direction of the fitted velocity gradient. The full version of this figure is available as Supporting Information to the online 
version of this article. 
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Table A3. Rotational properties of SCUBA clumps located with CLFIND. These properties have been found by non-linear least squares fitting of a linear 
velocity gradient to the first moment of the C^^O / = 3 ^ 2 data. Only significant detections are shown with cr^/^ ^ 3. The full version of this table is 
available as Supporting Information to the online version of this article. 



Sub-region ^ 


Clump ^ 


Hatchell ' 


Number ^ 


^±Gc^' 


%/ 


^/Gc^S 


iSrot^ 


j' 




ID 


class 


of points 


(kms-ipc-1) 


(degEofN) 






(km s ^ pc) 


NGC1333 


2 


I 


1117 


3.16±0.76 


155.44 


4 


6.98E-04 


1.05E-03 


NGC1333 


3 





2106 


4.85 ±0.35 


26.88 


14 


5.39E-03 


3.51E-03 


NGC1333 


5 


I 


2229 


1.54±0.41 


21.10 


4 


8.91E-04 


1.63E-03 


NGC1333 


6 





978 


4.66 ±0.84 


125.79 


6 


3.35E-03 


1.44E-03 


NGC1333 


8 


I 


330 


10.63 ±2.67 


-58.59 


4 


4.67E-03 


9.83E-04 



^ Name of sub-region map. 

^ Clump identification from Curtis & Richer (2010). 

^ H07 clump class, S=starless, 0=Class protostar and I=Class I protostar. ^ Number of points fitted in the map. 

^ Fitted velocity gradient and its error. 

f Angle of the fitted velocity gradient, east of north. 

^ Level of significance of the fit, a^/^. 

^ Ratio of the kinetic to gravitational energy. 

' Specific angular momentum. 

Table A4. Rotational properties of SCUBA clumps located with GAUSSCLUMPS. These properties have been found by non-linear least squares fitting of a 
linear velocity gradient to the first moment of the C^^O / = 3 ^ 2 data. Only significant detections are shown with a^/^ ^ 3. The full version of this table is 
available as Supporting Information to the online version of this article. 



Sub-region ^ 



Clump ^ 
ID 



Hatchell ' 
class 



Number^ ^±a^^ %^ 

of points (km s ~ ^ pc ~ ^ ) (deg E of N) 



/(T^^ 



fto, 



J' 

(kms~^ pc) 



NGC1333 


1 





800 


5.33±1.33 


135.67 


4 


6.90E-04 


7.86E-04 


NGC1333 


4 


I 


922 


4.93±1.16 


-79.61 


4 


1.22E-03 


7.41E-04 


NGC1333 


7 


I 


1222 


6.73 ±0.84 


-77.95 


8 


4.27E-03 


1.56E-03 


NGC1333 


8 





812 


4.25±1.11 


26.44 


4 


1.07E-03 


5.54E-04 


NGC1333 


9 





1374 


5.37 ±0.65 


-38.35 


8 


3.70E-03 


1.63E-03 



^ Name of sub-region map. 

^ Clump identification from Curtis & Richer (2010). 

^ H07 clump class, S=starless, 0=Class protostar and I=Class I protostar. ^ Number of points fitted in the map. 

^ Fitted velocity gradient and its error. 

f Angle of the fitted velocity gradient, east of north. 

^ Level of significance of the fit, cr^/^. 

^ Ratio of the kinetic to gravitational energy. 

' Specific angular momentum. 
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NGC 1333 1 



NGC 1333 4 



NGC 1333 7 



Figure A2. SCUBA 850 |xm clumps identified with GAUSSCLUMPS containing significant velocity gradients, cr^/^ ^ 3 in NGC1333, IC348 and L1448. 
The colour-scale denotes theC^^O/ = 3^2 line centre velocity, linearly scaled from the clump minimum (blue/black) to maximum (red/white). Contours 
of 850 jim flux density are overlaid (black) at 100, 200, 400, 800, 1600 and 3200 mJy beam" ^. A white/black circle marks the position of the clump peak 
emission with a white line positioned to show the direction of the fitted velocity gradient. The full version of this figure is available as Supporting Information 
to the online version of this article. 



© 2010 RAS, MNRAS 000,[l]{T7] 



